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We present a detailed discussion of the solution of Nambu-Gor'kov equations in superfluid nuclei, 
which provide a consistent framework to deal with the interplay between particle-hole and particle- 
particle channel, including the effects of the fragmentation of the quasiparticle strength and of 
the pairing interaction induced by the exchange of collective vibrations. The coupling between 
quasiparticle and vibrations is determined from the experimental polarizability of the low-lying 
collective surface vibrations. This coupling is used to renormalize the properties of quasiparticles 
obtained from a BCS calculation using the bare nucleon-nucleon interaction. We apply the formalism 
to the case of the nucleus 120 Sn, showing results for the low-energy spectrum and the quasiparticle 
i— l strength distribution in neighbouring odd nuclei and for the neutron pairing gap. 
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I. INTRODUCTION 

The key role played by pairing correlations in atomic nuclei was recognized soon after the development of BCS 
theory. The number of Cooper pairs involved is small, and one can study nuclear superfluidity within BCS in terms 
of specific orbitals lying close to the Fermi energy, interacting through an attractive nucleon-nucleon force. However, 
the properties of the particles are considerably renormalized by the coupling to the collective vibration of the system, 
which leads to an enhancement of the effective mass and of the level density close to the Fermi energy, and to a 
fragmentation of the particle strength observed in one-nucleon transfer reactions. While these phenomena have been 
extensively studied in normal nuclei [2(2], their consequences on pairing correlations have been much less investigated. 
Obtaining an accurate description of these effects requires to go beyond the BCS framework. A convenient formalism 
that allows one to consider the interplay between the particle-hole and the particle-particle channel is offered by the 
Nambu-Gor'kov equations, used in condensed matter to deal with strong coupling superconductivity [3]-[5]. These 
equations imply the calculation of energy-dependent normal and abnormal self-energies, take into account the pairing 
interaction induced by the exchange of collective vibrations between members of the Cooper pairs and lead to theo- 
retical predictions concerning the low-energy part of the nuclear spectrum; they also produce the structure elements 
needed for a consistent calculation of one- and two-nucleon transfer reactions [5] . 

In this paper we extend our previous investigations concerning the microscopic origin of pairing (cf. in particular 
[7j- |T2]). In Section In] we provide a detailed account of the formalism, comparing two different schemes for the 
dealing with the Nambu-Gor'kov equations, either solving first the BCS equations with the bare nucleon-nucleon 



force and then adding renormalization effects (prior scheme, cf. Section II A I , or dealing at the same time with both 



sources of pairing (post scheme, cf. Section II B ) . The coupling between quasiparticle and vibrations is computed 
according to the basic rules of Nuclear Field Theory (NFT) [T5H15] , Use is made of the collective model of Bohr 
and Mottelson to calculate the properties of vibrational states in the Quasiparticle Random Phase Approximation 
(QRPA) with a separable force, with a coupling constant chosen so as to reproduce the experimental properties of 
low-lying collective surface modes. In this way we avoid the use of free parameters in the calculation, which then 
only depends on the choice of the interaction adopted to produce the mean field. In most of the computations we 
shall use the SLy4 effective interaction, which provides a good reproduction of the overall mean field properties, and 
leads to a reasonable level density around the Fermi energy, after including renormalization effects. We shall derive 
gap equations which generalize the usual BCS exp ression, lend themselves to useful approximations and allow one 
to make contact with other studies. In Section III we shall solve the Nambu-Gor'kov equations in the case of 120 Sn, 
comparing the theoretical low-lying spectrum in neighbouring odd nuclei with experimental data derived from one- 
neutron transfer reactions. We shall also compare the results obtained with a self-consistent iterative solution of the 
Nambu-Gor'kov equations with the quasiparticle approximation. It is possible to clearly separate the contributions 
to the pairing gap associated with the bare nucleon-nucleon force and with the renormalization effects. We find that 
the two contributions have comparable magnitudes, confirming previous studies. We shall also briefly consider the 
role played by spin modes, which provide a repulsive contribution to pairing in the So channel, whose magnitude is 
however very uncertain. Conclusions and a hindsight will be presented in Section \P7\ The sensitivity of our results to 
various elements of the calculations is discussed in the Appendix. 

II. THE FORMALISM 

The object of our study are pairing correlations in superfluid nuclei taking into account both the contribution of 
the bare nucleon-nucleon interaction, and the many-body effects associated with the coupling between particles and 
vibrations (cf. Fig. [I]). 

Our approach will be based on two basic assumptions: 

(A) We shall start from a mean field obtained with a Hartree-Fock (HF) calculation with an effective force. Our 
quantitative results will then depend to some extent on the choice of the force adopted to produce the mean field, 
because pairing correlations are very sensitive to the position of the single-particle levels. In principle, one could aim at 
determining the parameters of a new effective force, by comparing the results obtained including renormalization with 
experimental data. This would require an extensive investigation, considering several isotopic chains. In this work 
we shall limit ourselves to the study of renormalization effects in the single nucleus 120 Sn, adopting existing effective 
forces which lead to a reasonable description of the available data - although we are aware that the quantitative 
aspects of our findings might be improved by a more careful determination of the mean field. 

(B) We shall not consider explicitly the renormalization processes affecting the vibrations of the system [10l IT6l 
117] , The main diagrams renormalizing the energy of the phonons are shown in Fig. plbl) (self-energy correction 
renormalizing the energy of the particle-hole transitions) and Fig. |2[b2) (vertex correction renormalizing the particle- 
hole interaction), while the diagrams shown in Fig. [2jdl) and Fig. [2ld2) renormalize the transition strength. The 
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FIG. 1. Basic diagrams taken into account in the present study, which renormalize the normal and abnormal self-energies 
obtained in mean field calculations: polarization (a), correlation (b), induced pairing interaction (c) processes. 



explicit inclusion of renormalization effects on phonons, on par with those on single-particle levels, represents an 
ambitious program that has been attempted only in a few cases for closed-shell nuclei [T5H2"0] . We shall instead take 
the view that such renormalization effects, if taken properly into account, should lead to agreement with experiment. 
In particular, we shall take the transition densities which are at the basis of the interweaving of single-particle 
(quasiparticle) degrees of freedom and collective modes from the collective macroscopic Bohr-Mottelson model. We 
shall determine the Particle- Vibration Coupling (PVC) strength (cf. Fig. [2Fb2)) from a QRPA calculation tuned to 
reproduce the experimental properties of the low-lying phonons, using a mean field with an effective mass m* = m 
that reproduces the experimental single-particle level density (cf. Fig. |2jbl)) (see Section III for more details). 

Within this frozen-phonon approximation, as a possible alternative, one might calculate the PVC microscopically 
making use of the transition density of the RPA phonons obtained from a self-consistent calculation performed with 
the same effective force employed to obtain the HF mean field. Calculations of this kind have been performed several 
times in the past for non-superfluid nuclei, with various degree of approximations pQ; see [21] for a recent calculation 
aiming at a complete self-consistence. The main drawback in this approach lies in the fact that QRPA leads to 
a relatively poor reproduction of the experimental energy and of the transition strength of the low-lying collective 
vibrational states in semimagic nuclei (for recent QRPA calculations cf. e.g. [22], [23]), which provide the main 
contribution to renormalization effects. This is consistent with the fact that processes beyond QRPA can strongly 
renomalize the phonon properties. We shall instead take the view that renormalization effects, if taken properly into 
account beyond the frozen phonon approximation, should lead to agreement with experiment. As a consequence, 
we shall take the PVC from a QRPA calculation tuned to reproduce the experimental properties of the low-lying 
phonons, using a mean field with an effective mass m* — m that reproduces the experimental level density. In other 
words, we shall assume that vertex corrections are effectively included in our PVC. 

The actual implementation of assumptions (A) and (B) will be presented below in the next Sections. Other, less 
essential assumptions will be adopted in our calculation in order to reduce the computational complexity. First of 
all, we shall limit the investigation of renormalization processes to states lying close to the Fermi energy. We shall 
then take into account the fact that the bare interaction - even considering its soft-core or V\ ow k versions - couples 
the single-particle states lying close to the Fermi energy with states lying up to several hundreds of MeV, while the 
phonon- mediated pairing interaction acts between pair of states separated by a few MeV [TTJ [23]. This suggests 
the convenience of a separate treatment of the two interactions in the pairing problem, by performing a HF+BCS 
calculation prior to the calculation of many-body effects. We shall refer to this approach as the prior scheme. 
Otherwise, one is forced to consider the bare and the induced interaction on the same footing (post scheme). In the 
following we shall first illustrate in detail the prior scheme, which will be used in most calculations presented in this 
paper. We shall then outline the main modifications involved in the post scheme. 







FIG. 2. The self-energy (a) of the vibrations calculated in the QRPA approximation is renormalized by the basic diagrams 
(bl) (leading to a renormalization of the particle-hole energies) and (b2) (vertex correction renormalizing the particle-hole 
interaction). Similarly, the QRPA transition strength (c) is renormalized by the corresponding diagrams (dl) and (d2). 



A. The prior scheme 



A convenient formalism for the calculation of the properties of quasiparticles in superfiuid nuclei within the prior 
scheme was given by Van der Sluys et al. [25 , although they did not devote particular attention to the renormalization 
of pairing correlations. In this approach one first accounts for the action of the bare force with a HF+BCS calculation, 
discussed in more detail in Section |IH A 1| , leading to quasiparticle energies E a , quasiparticles amplitudes u a and v a 
and to a pairing gap A^ cs . One then renormalizes the obtained quasiparticles including the coupling to vibrations 
calculated in the QRPA. The derivation of the formalism was based on the equation of motion method, which has a 
very close relation to the Green's functions formalism used to derive the Nambu-Gor'kov equations, commonly used 
to study superconductivity in condensed matter physics [26] . 

In order to calculate the renormalization of a quasiparticle in spherical nuclei, denoted by its associated quantum 
numbers a = {nlj}, one has to solve a system of linear equations obtained coupling the quasiparticle with more 
complex configurations including phonon states through the basic vertices V and W shown in Fig. [T] The phonons 
are characterized by their angular mom entum A and by their energy fiuj\^. We shall assume that phonons have natural 
parity, it — {— 1) A (see however Section III A 3 1. For illustration, we write the diagonalization problem including only 
two other quasiparticle states b and c, and a single phonon A, v. 
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Many eigenvalues -E a (n) and eigenstates a(n) are obtained from the diagonalization of the matrix lib, giving rise to a 
fragmentation of the associated quasiparticle strength. For a given eigenvalue E a ( n \ > there exists a corresponding 
eigenvalue E a t- n ) = —E a / n \. As in standard BCS theory, we will keep only positive energy solutions, n > 0. 
The amplitudes obey the normalization condition 
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C a (n),b,\v and D a / n \ b>Xl/ being the components on the complex lgp ® lph states, while cc r n ) and y Q („) are the 
components on the original lqp state a. The fragment a(n) carries a fraction of the strength 



N a(n) = X 2 a{n) + V 2 a(n) < 1 



(3) 



which is to be compared with the experimental quasiparticle strength, as determined for example in one-particle- 
transfer reactions. 

The excitation operator appearing in Eq.(l) 

°a(n) = a a{n) qp + C a(n)b\uOC b T Xv - D a{n)bXly a b T X p 

contains the quasiparticle component 

~t _ t 

a a{n) qp — X a(n)C*a Va{n)Ot a - 

Taking into account the BCS Bogoliubov transformation, 



Oil = U a a \ + v aO-a 



one can write 



where the quantities 
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V a (n) = Xa(n)Va + Va(n)U a - 



(4) 



represent the new quasiparticle amplitudes associated with a given fragment a(n): their squares give the spectroscopic 
factors associated with one-nucleon transfer reactions. 



We can then calculate the renormalized pairing gap as (cf. below Eq. (38)) 
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The quasiparticle-phonon matrix elements V(abXv) and W{ab\v) in Eq. (fTl) are given by 

V{abXv) = [|££] V2 E c < d (l + <U~ 1/2 x 
[X cd {\v)V{cd\b- a) + {-iy«-^+ x Y cd (\is)V(cd\a; b)] (6) 

and 

W{abXu) = [f±^] V2 £ c < d (l + ^)- 1/2 x 
[X cd (A^)i?(cdA6; a) + (-l)^-^+ A y cd (A^)Q(«fAa; b)] , (7) 

where A c <2 and Y cd are the forward and backward amplitudes resulting from the QRPA calculation. 
The terms V(cdXb;a), Q(cdXb;a) and R(cdXb;a) are given by 

V(cdXb;a) = -(u a v b u c u d - v a u b v c v d )G(abcdX) 
+ (u a u b u c v d - v a v b v c u d )F(abcdX) 
+ (-l) jc ~ jd+x (u a u b v c u d - v a v b u c v d )F(abdcX) (8) 



Q(cdXb; a) = {u a u b u c u d + v a v b v c v d )G(abcdX) 
+(u a v b u c v d + v a u b v c v d )F(abcdX) 
+ {-iy^ ]d+x {u a v b v c u d + v a u b v c u d )F(abdcX) (9) 



R{cdXb] a) — —(u a u b v c v d + v a v b u c u d )G{abcdX) 
+{u a v b v c u d + v a u b u c v d )F(abcdX) 
+ (-iy--^+ x ( UaVbUc u d + v a u b v c v d )F(abdcX), (10) 

where F(abcdX) and G(abcdX) denote the angular momentum coupled antisymmetrized particle-hole and particle- 
particle < a&A|V^|cdA > as matrix element respectively. Note that in the limit of non-superfluid nuclei (all u,v terms 
equal to or 1), the matrix elements V(abXv) connect pairs of states a, b above or below the Fermi energy through the 
F— terms, and pairs of states on opposite parts of the Fermi energy through the G— terms; while the opposite is true 
for the W(abXv) matrix elements. In the following we shall not take into account the coupling with pair vibration 
modes, neglecting the G terms in Eqs. ((8])-(10). While this coupling is known to be relevant for closed shell nuclei, 
it is expected to be much less important for the superfluid case, since most of the two-particle transfer strength is 
already incorporated in the BCS ground state (gauge space deformed) wave function J27) . 

The QRPA amplitudes could be obtained from a calculation performed using the same force adopted in the HF+BCS 
calculations. In our approach, the QRPA calculation is instead decoupled from the renormalization process, in keeping 
with the main assumption (B) discussed above. In our QRPA calculation we shall use the separable force 



BTJ r)JJ 
V{r 1 ,r 2 ) = -K sel} r 1 —r 2 —Y J XxY^{6 l )Y x ,{6 2 ) (11) 

Xf^i 

where U(r) is a potential that gives a good reproduction of the experimental levels. In practice, we adopt the 
Woods-Saxon parametrization given in [28_ (cf. Eq. (2-182)) together with an empirical pairing coupling constant 
adjusted to reproduce the pairing gap deduced from the experimental odd-even mass difference. The parameters xa 
are determined so as to get a good agreement with the observed properties (energy and transition strength) of the 
low-lying surface modes. More precisely, we shall reproduce the polarizability f3 xi /ft,uj\i of the low- lying modes, where 
fi\ v denotes the experimental nuclear deformation parameter. In fact, the matrix elements of the phonon- induced 
pairing interaction for levels close to the Fermi energy are approximately proportional to the polarizability of the 
mode (cf. Eq. ( |30| below). The resulting values of xx turn out to be close to 1 (cf. Section [TtT| ) , indicating that the 

QRPA coupling constant is close to the Bohr-Mottelson self-consistent coupling constant K se if = / r^r^r 2 dr 

This scheme then reduces to the collective particle-shape vibration (phonon) coupling scheme given by Bohr and 
Mottelson 15J (cf. Eqs. 6-207- 6-209). In fact the particle-hole matrix elements, neglecting the exchange terms (cf. 
on this point [53], Eq. (14.54) and Chap. 16), are given by 

F(abcdX) = -KseifXx < abXfi\ri^-Y^(6i)\0 >< Q\r 2 ^-Y^{6 2 )\cdX^ > (12) 



where /j, is any of the z-projections of the angular momentum A. In this expression the QRPA-like single-particle indices 
(c, d) and the scattered particle indices (a, b) appear in separated factors, so that one gets the angular momentum 
reordering property F(abdcX) = (-lp<=-^+ A F(a6cdA) = (-l) j "- jb+x F(bacd\), and 
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The quantity in the summation is precisely the transition amplitude M{Xv) of the M = r 2 ^rY x (62) operator, which 
is usually expressed in terms of the so-called collective deformation parameter as M(Xv) — ol° Xv j K se if, assuming a 



collectively deformed density dp 
In this way we can write 
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Finally, following the notation in [TS], Eqs. (6-207 to 6-209) using the reduced matrix element < j&| |^a| |j a >= 
(— l) Ja ~i h < j a jb] \fj,\Y\p\0 > \/2A + 1 and the relation a° Xl , = (i\ v /\/2\ + 1, we can write 



V(ab\v) — h(abXv)(u a Ub — v a Vb), 



(15) 



where 



which is the basic vertex in 
Analogously one finds 
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corrected by our effective deformation parameter /9t = X\vfi\v 



W{ab\v) = h{ab\v)(u a Vb + v a Ub) 



(16) 



(17) 



A self-consistent renormalization procedure by iterating the diagonalization process, using the previously renormal- 
ized quasiparticles to build the complex lqp ® lph states. This means that the basic V, W matrix elements are now 
calculated as 



V{ab{m)\v) = h{ab\v)(u a u b ( m ) - v a v b ( m )) 
W(ab{m)\u) = h{ab\v){u a v b { m ) + v a u b ( m )) 



(18) 



where the u a and v a amplitudes obtained from the initial HF+BCS calculation are kept fixed in the iteration process; 
on the other hand the amplitudes Ub(m) an d Vb( m ) refer to the amplitudes associated with the to— fragment resulting 
from the renormalization of the state b in the previous iteration step. In this way a consistent fragmentation of the 
different states is constructed through the iterative procedure. 

The iteration process gives rise to the so called no-line crossing rainbow series of the self-energy (see below), in 
which all orders are summed up coherently, and which is expected to play an important role in the limit of strong 
coupling (cf. FigJ3j). The size of the matrix to be diagonalized increases exponentially at each iteration, and some 
sort of numerical approximation is needed, as will be discussed in Section |HI| 



1. Energy dependent, BCS-like formulation 

The diagonalization of the eigenvalue problem (nj is equivalent to the 2x2 energy-dependent eigenvalue problem 
derived from an approach based on Green's function formalism (cf. e.g. [8]): 
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FIG. 3. The coupling of the quasiparticle to many-phonon states, corresponding to the rainbow diagrams shown in the figure, 
is included in our approach through the iteration of the self-consistent diagonalization of the Nambu-Gor'kov matrix. 



where one has introduced the energy-dependent, normal self-energies E Q A° and E , p -.° given by 
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and the abnormal self-energy 
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We note that E 11 evaluated at a given energy E is equal to — E 22 evaluated at — E, that is, S A° 
normalization of the quasiparticle strength of the n- fragment (cf. Eq. Q) is given by |30j 
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In this form one can easily make contact with the formalism based on the Green's function for superfiuid systems, 

G a (E a + id) = {(E a +iS)l- E a r 3 - tf°(E a + id)}' 1 , (23) 

where T3 denotes a Pauli matrix, and where 
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is our phonon-mediated self-energy matrix evaluated at the complex energy E a + iS, which coincides with the matrix 
introduced by Nambu and extensively used in condensed matter to deal with strong coupling superconductivity 26J. 
Eqs. ( 23 ) , |24|) must be solved self-consistently due to the fact that the quasiparticle strengths needed for the evaluation 



of Y.P ho (E a + iS) through the V and W matrix elements are obtained from the imaginary part of G a (E + iS) (cf. [25], 
Eqs. (46-47)). 

In order to get more insight concerning the respective contributions of the bare and of the phonon-induccd interaction 
to the pairing gap, it is useful to rewrite the 2x2 eigenvalue problem (19) in terms of the amplitudes u, v, instead of 
x, y, by inverting the relation Q, obtaining 
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where e a denotes the HF single-particle energy, while the new normal self-energies are given by 
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One can separate the abnormal self-energy into two terms, writing 
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The first term, A^ cs , is the pairing gap associated with the bare interaction obtained in the HF+BCS calculation, 
while the second term is associated with the phonon-induced interaction and is given by 
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Using Eqs. (18), (20) and (21) this expression can be simplified and the abnormal self-energy can be rewritten as 
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where we have introduced the induced pairing interaction: 
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Furthermore, we can symmetrize the matrix ( 25 ) in order to get a 2 x 2 eigenvalue equation which is formally identical 



to the BCS eigenvalue equation. This can be achieved multiplying Eq. (25) by the Z a i n ) energy-dependent function, 
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We note that according to the definition above, Z is the inverse of the correspondent quantity as defined in [8].|26j. 
We also note that the symbol Z is often used instead of N to define the quasiparticle strength. In fact the two 
quantities tend to take similar values close to the Fermi energy (cf. Fig. [7] below). This similarity can be explained, 
on the one hand, by noting that for the lowest pole, the finite difference in the expression of Z can be approximated 
by a derivative: 



Ja{n) 



a(n) a( — n) 



IE, 



a(n) 



dE a{n) 



(33) 



On the other hand, in the normalization of the strength of quasiparticles close to the Fermi energy, x » y so that, 
using Eq. (p2|, 



n, 



a(n) — x a (n) ' Va(n) ~ X ' a(n) 



DY 11 v 



8E, 



It is possible to rewrite Eq. ( 25 ) as 



e a (n) - e F 



A 



a(7i) 



A a ( n ) -(e a ( n ) - &f) J \ V a(n) 



*a(n) 



a(n) 



Ea(n) ( ~ Q(n) 



where 



£o(n) ~ e F - Z a{jl) (e a -e F ) + E^™ 



(34) 



(35) 



(36) 



and where Ejj/*? is the even part of SH > : 
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yll _ y*22 yllpho _ y22pho 

^even _ a (") a (") __ / ,2 „,2\ a (") a (") 

^(n) — 9 — \ U a ~ v a) 9 



2u„.v„Y l 



12pho 
a(n) 



(37) 



The term e a („) in Eq. (35) represents the renormalized single-particle energy, and one can now identify the pairing 



gap with the term A a („) |26| 



A a („) — Z a ( n )E a(n) — Z, 



a(n) \& a 



BCS , fil2,p/io\ _ Xb. 



a(n) 



A 



a(n) 



^pho 
a(n) 



(38) 



This quantity coincides with the expression[5]introduccd above, which can be used when solving the energy independent 
problem R\, because it involves only quasiparticle energies and amplitudes. Note also that the quasiparticle energy 
relates to the new gap and single-particle energy as in the usual BCS equations: 



E. 



'a(n) = sJ(Za(n) ~ &f) 2 + A 2 (n) 



(39) 



It is now possible to write a generalized gap equation. Introducing the total quasiparticle strength for a given 
fragment, 



the product u a (n)V a (n) rnay be obtained in the usual way from the 2x2 secular equation as 



(40) 



U a (n)V a (n) = N a ( n ) 



A 



a(7i) 



2j(e a{n) - e F r + A 2 o(n) 



(41) 



U a (n)V a (n) = N a ( n y 



Substituting in Eq. ( 29 ) one obtains 



7 , n't 12 



^(„)( e «-^+ S aw) 2 + (^(") S a(n)) 



(42) 



fil2 



A„ 



E 

6(rn) 



V ind {anbm)N b ( m ) 



V, 12 



2J(e b - e F 



n b(m) 



en\1 



(S 12 



b(m 



) 2 



(43) 



The second term in this equation is a generalization of the usual BCS gap equation, and clearly demonstrates how 
the action of the different fragments of the original quasiparticles is modulated by their quasiparticle strength N b i m \ . 
The equation, however, is of little practical use as it stands because it involves the energy dependent interaction Vi n d 
which contains a "dangerous" denominator (cf. Eq. (30)). The formula will be further discussed in Section III A 4 
and we shall present a similar expression in Section II B (cf. Eqs. 53p4 ). 

The approach presented so far (and in [25) ) neglects the lgp-exchange interaction between the complex lqp® lph 
states (see fig. 6.10 c-d in |15j). In fact the associated matrix elements are set to in the matrix (fTl). The ignored 
processes would account for some violation of the Pauli principle arising from the microscopic structure of the QRPA 
phonons, which may imply double occupation of the quasiparticle state in the com plex l qp ® lph state. These 
violations are anyway very small for the calculation reported in this paper (cf. Section VI D ) . 



Those processes account also for vertex renormalization terms in the self-energy, which are not taken into account 
in the rainbow series we have considered. While in condensed matter they are usually neglected based on the Migdal 
theorem, in nuclear physics they are usually considered to be of minor importance because their contribution to the 
self-energy implies a recoupling of angular momenta that when summed up over all possible intermediate states is 
expected to lead to a rather strong cancellation [3TJ[35]. As we have mentioned discussing approximation (B) above, 
we shall assume that they are implicitcly included in our effective PVC, computed making use of phcnomcnological 
phonons and single-particle levels. 

Although within NFT tadpole diagrams should be included |33j . we will neglect the energy- independent contribu- 
tions associated with them (cf. Fig. Pfl) , which take into account the effect of zero-point fluctuations on the quasiparticle 
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FIG. 4. Lowest order diagrams representing the renormalization of the quasiparticle self-energy by zero-point fluctuations. 

energy [36] . This kind of diagrams modifies the nuclear density (cf. Fig. O and plays a particular important role in 
the calculation of nuclear radii [3TH3UJ , representing the leading correction beyond mean field for closed shell nuclei, 
and producing sizeable isotopic effects [3D] . However, they lead to relatively small changes in the mean-field potential. 
The shifts of the single-particle energies for A = 120 can be estimated to be of the order of 150 keV [41], and, being 
of static nature, we shall assume that they are effectively taken into account in the mean field. The tadpole diagrams 
can also influence the abnormal density an the calculation of the bare pairing gap |42) , but no calculations have been 
performed for superfluid nuclei. We have estimated that the effect of diagrams (a) and (b) in Fig. p3] changes the 
abnormal density k by about 5%. This value can be obtained from the expression 



Sn(r) 



ab\ 



(2X + l)Y^ bX K a (r)/(2j a + l), 



(44) 



where Y a b\ denotes a backward amplitude calculated in QRPA. The diagrams (c) and (d), in the normal case, give 
a contribution to the renormalization of the radius which is about three times larger than the one produced by (a) 
and (b) J37J. This can be considered as an upper limit since the strong r 2 dependence tend to enhance their relative 
importance. 

Thus we estimate a total tadpole effect to the pairing field of less than 20/value, that is less than 200 keV, which 
is a number rather consistent with that estimated above for the normal single-particle energy. 

Finally, we note that the formalism presented in this Section is similar to the one adopted by Avdeenkov and 
Kamerdzhiev 34, 35 . However, a precise comparison of their results with ours is difficult, because they have followed 
a different approach, trying to extract 'bare' single-particle levels and pairing gaps from the experimental levels, 
instead of renormalizing the levels obtained from an effective mean field. 



B. The post scheme 



In the previous Section we have outlined the prior scheme in the calculation of the pairing properties of superfluid 
nuclei. That is, we have started from a HF+BCS calculation which accounts for the bare pairing interaction and 
we have then added the effects of the PVC diagonalizing the matrix in Eq.(fi]). However one could also adopt a 
different scheme (which we will call the post-scheme) in which one starts from the HF solution and then one includes 
simultaneously both renormalization effects and the bare pairing interaction through the iterative procedure. This 
approach is particularly appropriate when the phonon induced pairing provides the leading contribution to pairing 
correlations, and is certainly needed when the bare interaction alone is not able to produce a superfluid solution. This 
is the case, for instance, in calculations of halo nuclei based on the PVC coupling [351133]. On the other hand, the 
prior-scheme is more natural, in the framework of general schemes based on the corrections to mean field properties. 

In the following we shall apply the post scheme to the calculation of pairing properties in 120 Sn, which is a well 
bound, superfluid nucleus. In this case we expect that the two schemes should give similar results concerning physically 
relevant quantities, namely quasiparticle energies, spectroscopic factors and pairing gaps. 
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FIG. 5. Lowest order diagrams representing the renormalization of the density operator (represented by the cross) by zero-point 
fluctuations in closed shell nuclei. 

The calculation in the post scheme is performed by introducing the bare pairing field £ 12bare associated with the 
bare interaction in Eq.pJ), by writing 



/ E a V{ab\v) V{ac\v) W(ab\v) W{ac\v) 

V{ab\v) huj X u + E b 

V{ac\v) hujxv + E c 

W(abXiy) -Uu\ v - E b 

W{ac\v) -huxu - E c 

\=FS 126are W{ab\v) W{ac\v) -V(abXv) -V{acXv) 



-r-y^l2.bare \ 

W(ab\u) 
W(acXv) 
-V(abXv 

—V{acXv 
-E a J 



I X <n) \ 



a 



c, 



a(n),6,Ay 



a(n),c,Ay 



-D 



'a(ri).b,\v 



-Va(r. 



( X <n) \ 



E, 



a(7i) 



c, 



a(n),c,Xu 



a(n),6,Af 
a(n),c, \u 



\ -Va{n) I 



(45) 



where now the unperturbed " quasiparticle" energy E a contains no contribution from pairing and is simply equal to 
the difference between the HF single-particle energy and the Fermi energy, 



E, 



a c -a 



£F , 



where S obeys the equation 



Z 12bare = -J2V bare (ab) 



(2J6 + 1). - 
S Ub(n)Vb(n) 



(46) 



b , n 
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and the -(+) sign is to be used for particle- (hole-) states, that is e a > ep (e a < ejr). 

The amplitudes u (n) i w a (n) are given by Eq. (H|, taking into account the fact that the initial u a and v a factors to 
be inserted in the iterative solution of Eq.Q are now equal to 1 or 0, depending on whether the level is above or 
below the Fermi energy: 



^a(n) •^a(n) ^a(n) Va(ji) 



(47) 



and 



Jj a(n) 



Va{n) ^a(n) %a(n) 



(48) 



respectively in the case of a particle- and a hole-state. 



Note that, putting the V — W matrix elements to zero, Eq. (46) reduces to the standard BCS equation for the 



pairing gap. The solution must be obtained through a self-consistent iterative procedure which in the general case 
involves simultaneously the bare interaction and the PVC vertices V and W. In general this is not equivalent to the 
two-step diagonalization performed in the prior scheme, in which one considers first only Y} 2bare (which is exactly 
equivalent to solve the bare BCS problem), and then takes into account the matrix elements T^'s and Ws. Nevertheless 
we expect that the two schemes lead to similar results, as long as the first diagonalization leads to a finite value of 
the order parameter, the pairing gap. 

As in the prior scheme one can use the 2x2 energy dependent BCS-like matrix to solve the problem. The only 
difference with respect to Eq. (35) is that now the renormalized pairing gap must include a term 5] 12bare representing 



the contribution of the bare pairing interaction: 



A„ 



ry l v-\126are 



\12,pho 



This expression can be rewritten in the more appealing way 

(2J6 + 1) 



x a(n) 



' pho 
a(n) ' 



A 



a(n) 



-z a(n) y^ 



V e ff(a(n)b(m))ub(m)Vb(r 



(49) 



(50) 



where one has introduced the effective interaction V e //, which is the sum of the bare and phonon- induced interactions: 

V e ff(a(n)b(m)) = V bare (ab) + V ind (a(n)b{m)) (51) 

Furthermore one can obtain a very closed form for the gap equation by eliminating the amplitudes u and v using 
Eq. p5b: 



U a (n)V a (n) = -Wo(n) 



»a(«) 



IE, 



N, 



a(n) 



a(n) ' 



7 , ,V 12 
Z 'a(") Zj a(n) 



2JZ 



a(n)( £ a ~ e F + ^^) 2 + (Z a(m) S^) 2 



(52) 



leading to 



A („) - Z a(n) E a(n) - 



E b(m) ? ^V;//(a(n)6(m))iV 6(m): 



J a(n) Z^b(m) 



y/(<*-eF+Vr l %) a + &% m) ) 2 ' 



(53) 



We note that for levels near to the Fermi energy the Z— factor is very close to the quasiparticle strength N. A 
more compact expression, bearing a direct resemblance to the standard BCS gap equation, may be obtained by 
reintroducing the Z— function both in the numerator and in the denominator: 



M") 



-Z a (n) Yl y e//(a(n)6(m))iV, 



A. 



b(m) 



b(m) 



b(m) 



2E, 



(54) 



b(rn) 



This equation is an exact consequence of the Nambu-Gor'kov energy-dependent problem and can be used as an useful 
starting point for approximate gap equations (see Section [ill A 4|) . In particular, restricting the sum to the main peak 
rn = 1 and neglecting the difference between N and Z, Eq. ( |54[ ) formally reduces to an approximate gap equation 
previously presented in the case of uniform matter ([3SJ (cf. Eq7(12)), [IS]). 
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Exp. 


Theory 


A" 


hujxi 


flu 


/3 2 xl /hu) X1 


Tiuxi 


^Al 


(Pt{ f f/hcoxi 


XA 


2+ 


1.17 


0.13 


0.014 


1.20 


0.13 


0.013 


0.86 


3" 


2.40 


0.16 


0.011 


2.71 


0.16 


0.010 


0.95 


4+ 


3.10 


0.11 


0.004 


2.33 


0.10 


0.004 


0.91 


5" 


2.27 


0.08 


0.003 


2.48 


0.09 


0.003 


0.91 



TABLE I. In the first columns we list the experimental energies htuxi in (MeV), deformation parameters fixi and polarizabilities 
Pxi/f l <^\i °f the low-lying states associated with the 2 + , 3~, 4 + and 5~ multipolarities. They are compared with the corre- 
sponding quantities calculated in QRPA, making use of the effective deformation parameter f3 x J — x>.Pxv In the last column 
we give the values of xa- In the case of 4 + , the experimental low-lying strength is fragmented in four peaks lying between 2.2 
and 3.8 MeV; the numbers in the table refer to an average weighted with the transition strength of each peak. 



III. 



RESULTS 



In the following we present the solution of the Nambu-Gor'kov equations for 120 Sn in various approximations. In all 
cases, we shall limit our attention to the five neutron levels belonging to the shell around the Fermi energy (namely, 
the l<77/2, the 2d 5 / 2 , the 3s 1 / 2 , the 2d 3 / 2 and the lh n / 2 orbitals). Most of the calculations will be performed in the 
prior scheme, making use of the Argonne nucleon-nucleon interaction as the bare pairing force in the ^o channel, 
which gives by far the dominant contribution to T = 1, J = pairing |47j . In Appendix VI B we shall also show a few 
results obtained with the Vi ow fe potential. 

In general we shall iterate the renormalization process. This requires looking for a self-consistent solution, either by 
successive diagonalizations of the matrix (11) or ( 45 ) , or looking for the solutions of the equivalent energy-dependent 
2x2 problem (19 1. Nearly all of the results presented in this work have been obtained from the the diagonalization 
procedure. However, we have verified in several instances the numerical agreement between the two approaches. In 
order to control the iteration process, one must introduce a cutoff procedure, or perform some averaging, in order 
to avoid the exponential increase of the number of solutions, retaining at the same time the essential information 
[3B],[3H]- We shall make use of two simple numerical procedures. A cutoff procedure can be adopted when one 
solves the energy-dependent problem. In this case one can limit the number of solutions kept at each iteration, by 
selecting only those fragments carrying a quasiparticle strength larger than a given cutoff N cut . An extreme case 
is represented by the so-called one-quasiparticle approximation, in which one keeps only the most important pole 
for each orbital a. We shall instead make use of the averaging procedure when we solve the energy-independent 
problem diagonalizing the matrices dTl or ( 45 1 . In this case we define a number N zones of energy zones. The zones are 
generally not of the same size, but are chosen so as to reflect the main features of the quasiparticle strength function. 
After each iteration we collect all the strength obtained within a given energy zone into a single fragment, placed 
at the calculated average energy position. In this way, the number of solutions associated with a given orbital a is 
kept fixed to Ndi m = 2x [N p i lon x N sps x N zones + 2] , where N p h on is the number of RPA solutions retained in the 
diagonalization and N sps is the number of single-particle levels considered. 



We shall compute the PVC coupling by first performing a QRPA calculation with the separable force (11) for 
the multipolarities and parities A 71 ' = 2 + ,3~,4 + and 5 _ . In each case, the coupling constant as determined so as 
to reproduce the experimental value of the ratio B(E\) /fiuj\^i (cf. the discussion about the approximation (B) in 
Section In] above). Experimental data are taken from [50]. The quasiparticle states used in the QRPA calculation are 
obtained from a BCS calculation with a monopole pairing interaction —GqP^P based on the levels of a Woods-Saxon 
potential parametrized as in ref. |28j . Eq. (2-182). The pairing coupling constant Go is adjusted so as to reproduce 
the value A exp ss 1.4 MeV derived from the experimental odd-even mass difference. The PVC matrix elements are 
then obtained from Eq. (16). 



Concerning the QRPA spectrum, we adopt an averaging procedure similar to that adopted for the quasiparticle 
strength. We include explicitly the strong collective low-lying vibrational states and we collect the remaining strength 
in a small number of peaks, reflecting the gross structure of the response for each multipolarity [5]. We have verified 
that the results are not sensitive to the details of the high-lying phonons. The properties of the low-lying phonons 
employed in the calculations are listed in Table [TJ where they are compared with the available experimental data. 
We remark that the values of xx (cf- Eq. (11)) are close to 0.9, reflecting the rather good accuracy of the collective 
model. 
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A. Calculations with bare pairing potentials 

In this Section, we shall discuss solutions of the Nambu-Gor'kov equations making use of bare nucleon-nucleon 
potentials as pairing interactions. As we have remarked above, we shall adopt the prior scheme. The practical 
difficulty with the post scheme is that the basis needed to account for the realistic bare interactions must span a 
broad energy interval (about 1 GeV in the case of the Argonne Vu interaction, and about 200 MeV in the case of 



Viowk)- This implies a large numerical effort in order to handle Eq.(45l, in particular because the matrix associated 
with a given angular momentum must include states with different number of nodes. We shall further consider the 
relation between the prior and post schemes presented in Section [XTJ making use of a simplified bare interaction. 

1. Bare pairing gap 

The first step in our calculation is represented by the solution of the HF+BCS equation with the bare force in the 
pairing channel; we shall not consider the influence of pairing on the mean field. 

We notice that the calculation with the bare force is in fact an extended BCS calculation, because for a hard-core 
interaction it is essential to consider the coupling between pairs of levels with different number of nodes, in order to 
take properly into account the coupling between bound and continuum levels [5TJ [52] . We take into account states 
up to 1 GeV. As a consequence, for a given set of quantum numbers a = {lj} one obtains a set of quasiparticles 
{Ek}', to each quasiparticle is associated an array of quasiparticle amplitudes u n k and v n k, which are the components 
of the quasiparticle states over the HF basis states 4> n ij- Going to the canonical basis, where the density matrix 
takes a diagonal form [53] . we look for the state having the largest value of the abnormal density, u rnax v max . As a 
rule, for a stable nucleus like 120 Sn, this canonical state is dominated by the quasiparticle state having the lowest 
value of the quasiparticle energy, E m i n . We then approximate the extended BCS calculation associating the value 
of v max and E m i n , treating them as proper BCS quantities; in particular we derive an associated pairing gap as 
A = 2u max v max E m i n , which is very close to the diagonal value of the gap in the original basis. The value of 

Vmax ? E min and A BCS are then employed as input values for the solution of the Nambu Gor'kov equations in the prior 
scheme described in Section I1A (cf. Eqs. (fTl) and Q, where they are denoted u a ,v a and E a ). This approximation 



leads to a substantial simplification of our numerical scheme. For nuclei close to the line of stability, where the 
canonical states have significant weight on continuum state and the PVC easily connects bound and continuum 
states, this approximation would not be justified and one should rather consider different numerical schemes, for 
example based on the continuum Green's function. 

The mean field will be taken from a HF calculation with the SLy4 interaction [54]. This interaction gives a good 
reproduction of the bulk properties of nuclei; moreover the resulting level density close to the Fermi energy (associated 
with an effective mass m^ rj 0.7m), once increased by renormalization effects, is in reasonable agreement with the 
experimental one (the associated effective mass increasing to to* = to^to^ rj to). The energies of the five single- 



particle levels lying closest to the Fermi energies are shown in Appendix VIA (cf. Fig. [19]). The detailed features of 
our calculations are of course influenced by the specific properties of the mean field, and in particular by its effective 
mass. In order to have some insight about the dependence on the adopted mean field, in the Appendix |VI A| we 
provide results obtained making use of two different Skyrme interactions, having effective masses smaller and larger 
than the SLy4 interaction. One should distinguish two different effects that influence the pairing gap calculated with 
the bare pairing interaction: on the one hand, the pairing gap is sensitive to the detailed position of the levels close to 
the Fermi energy, that could be influenced by static contributions not considered here, as those produced from tensor 



correlations J55j[56] and from tadpole diagrams [42] (for the latter, cf. the comments at the end of Section II A); on 
the other hand, the pairing gap is also sensitive to the momentum dependence of m^ for large momenta. 

In a previous work [10J we presented a solution of the HFB equations using the Argonne vn equation in the pairing 
channel, which led to a pairing gap of the order of 700 keV. There we used a modified SLy4 mean field, reducing the 
spin-orbit coupling strength by about 15% (furthermore we included in the Hamiltonian the terms proportional to 
the square of the spin density, which will not be included in the following, and used a slightly larger PVC strength). 
This was done, in order to improve the agreement of the results obtained after renormalization with the experimental 
spectrum and pairing gap. Here we prefer to take the HF field obtained with the original SLy4 interaction. This is 
associated with a larger level density around the Fermi energy. As a consequence, the obtained pairing gap, shown in 
Fig. [6] (full dots), is larger, being equal to about 1.1 MeV. Similar calculations have been performed with the V\ ow k 
potential using the same mean field [ST]. In that case, the result depends on the cutoff A adopted to obtain V\ OVJ k'- 



for values of A smaller than about 4 fin. , the gap is close to 1.4 MeV, (cf. Fig. 25), while for increasing values of A 
the Viowfe potential reduces to the Argonne potential and one obtains the result shown in Fig. [6] We notice that our 
numerical results with the Argonne potential are in very good agreement with those reported in ref. [57] for very large 
values of A. The reason for the difference between the results obtained with the Argonne and V\ ow k interaction has 
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FIG. 6. The state-dependent neutron paring gap A calculated in BCS with the bare v%4 interaction is shown as a function 



of the SLy4 HF single-particle energy of the five valence orbitals (cf. Fig. 191, and is compared to the renormalized gap A 
(cf. Eq. ( 38 1 ) obtained solving the Nambu-Gor'kov equations by iteration; also shown (dashed line) is the gap obtained at 
the first iteration. We also show the decomposition of A into the bare and phonon contributions A bare and A pho . The value 
of the Fermi energy ep and of the gap obtained from the experimental odd-even mass difference /S. exp are also indicated. 



been discussed in detail [57. 58 : using the Argonne potential implies using the Skyrme effective mass mu (equal to 
about 0.7m inside the nucleus) up to very high momenta, while the renormalization process leading to V\ ow k implies 
that mk goes to the free mass for momenta larger than the cutoff A. The lower value of rrik at high momenta leads 
to a smaller level density and to a reduction of about 300 keV in the value of the gap with the Argonne interaction 
with respect to the values obtained with V\ ow k for cutoffs A up to about 4fm~ 1 . An analogous difference shows up in 
calculation of the gap in uniform matter. While the lack of momentum dependence of the effective mass associated 
with Skyrme interactions is certainly unrealistic, the proper momentum dependence is still to be established; and the 
density dependence of rrik associated with the SLy4 interaction in nuclear matter is not far from that resulting at 
the Fermi energy from calculations based on Bruckner theory |58j . Furthermore, a precise determination of the value 
of the bare gap within the Viowfc scheme requires the consideration of the effects of the three-body force, which is 
expected to provide a repulsive contribution [59] . A calculation of three-body effects within the Viowfc renormalization 
scheme [60], [61] leads in fact to a reduction of the gap, down to values close to those obtained with the Argonne 
interaction in Fig. [6j An average value of the bare gap close to 1 MeV was also derived in the analysis of refs. [34JI35J. 



We have brought circumstantial evidence which testifies to the fact that a pairing gap A BCS sa 1.1 MeV for the 
levels around the Fermi energy obtained using the Vu potential as bare pairing interaction and the effective mass 
associated with SLy4 represents a reasonable starting point, being well aware that the determination of the mean 
field and of the associated effective mass is one of the most important issues which remains to be fully clarified, for a 
quantitative and microscopic calculation of the gap in finite nuclei. In Section [VI A| we investigate the dependence of 
the results on the adopted mean field, whille in Section [VI B| we provide some results obtained adopting V[ ow k as bare 
pairing interaction. 
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2. Solution of the Nambu-Gor'kov equations 



The state-dependent pairing gap obtained from the solution of the Nambu-Gor'kov equations is compared to the 
bare gap in Fig. |6j Renormalization effects lead to a total gap A about 600 keV larger than A BCS . Most of the 
effect is obtained already with the first diagonalization of the Nambu-Gor'kov matrix. The self-consistent iteration 
process leads to a further increase of the gap by about 10%. We recall that one can distinguish two contributions to 
the renormalized gap A = Z[A BCS + t, 12 ' pho ] = A bare + A pho , associated with the bare and with the phonon-induccd 
interaction (cf. Eq. p8|). They are also shown in Fig. [6] and turn out to be of similar magnitude. This confirms 
that the phonon-induced pairing interaction is as important as the bare interaction in determining pairing properties 
of heavy nuclei. Notice that a proper comparison between the role of these two sources of pairing can only be made 
analyzing their contribution to the final physical result, and not just by comparing the BCS and complete pairing 
gaps. In fact, processes beyond mean field act in a complex way, not only giving rise to the induced pairing interaction, 
but also reducing the effect of the bare interaction through the Z— factors. The values of Z for the five orbitals are 
shown in Fig. [7] where it can be seen that they are close to 0.7, bringing the bare contribution from the BCS result 
(about 1.1 MeV) down to about 0.8 MeV, which is about one half of the total renormalized gap. The other half is 
provided by the pairing induced interaction. The values of Z are similar to the values of the quasiparticle strength 
N = U 2 + V 2 , except for the orbital ^5/2- The values of N are in rather good agreement with the quasiparticle 
strength measured in one-nucleon transfer reactions, shown by stars in Fig. [7](cf. also below Figs. [9]and 10 with the 
related discussion), although one has to consider that the experimental values are affected by a large error, estimated 
to be about 30%. 

The renormalized pairing gap exceeds the experimental value obtained from the odd-even mass difference by about 
300 keV, the value of A exp lying in between A BCS and A. We remark that the estimate of A exp varies by about ± 
100 keV depending on which odd-even mass formula is used [52J H3] ( the five-point formula A' 5 ) yielding a value 
of 1.39 MeV), and that a more consistent comparison would imply a theoretical calculation of the binding energies. 
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FIG. 7. Comparison of the N~ and Z— factors associated with the lowest quasiparticle peaks in the Nambu-Gor'kov calculation 
shown in Fig. [m Also shown by stars are the values of the experimental quasiparticle strength [66], except for the d 5 / 2 orbital 
which shows a pronounced fragmentation (cf. text). 
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FIG. 8. The theoretical quasiparticle spectra obtained at the various steps of the calculation are compared to the experimental 
data. See the text for more explanation. 



Such a precise comparison of the gaps is probably not very significant at the present stage of the theory, given the 
uncertainties which affect both the BCS calculation - like the dependence on the adopted HF mean field and the 
effect of three-body forces mentioned above - and the renormalization process - in particular the exchange of spin 
modes, not included in this calculation (see the discussion in Section III A3 and Section IV A). However, the fact that 
the value of the renormalizcd pairing gap is found to be relatively close to experiment is satisfactory. The necessity 
of going beyond mean field turns out anyway to be clearer considering other physical quantities, in particular the 
energy spectrum of neighbouring odd nuclei with the associated strength functions and spectroscopic factors. This 
can be seen in Fig. [8] where we compare the features of the quasiparticle spectrum obtained at the different steps 
of the calculation. In the first column (HF), we report the absolute value of the difference \e a — ep\ obtained in the 
HF calculation with the SLy4 force, referred to the value of this difference for the level d 3 / 2 , which lies closest to the 



Fermi energy €f = -8.05 MeV (cf. Fig. 19). In the second column (BCS), we give the values of the quasiparticle 
energies obtained in the HF+BCS calculation with the Argonne force, referred to the lowest quasiparticle. In the 
third column (NFTO) we show the spectrum obtained taking into account processes beyond mean field calculated 
with the Nambu-Gor'kov equations without iterating, while in the fourth (NFT) we show the self-consistent solution: 
in these cases we show the energy of the lowest fragment for each quantum number (this is also the one collecting the 
largest quasiparticle strength, except for the d 5 / 2 orbital, which is very fragmented and whose case is discussed below): 
Finally in the fifth and sixth columns we give the position of the lowest peaks measured in 119 Sn and 121 Sn. The main 
discrepancy of the renormalized spectrum concerns the 7/2 + state, whose position is brought too low by the PVC by 
about 300 keV ( this is probably related to the initial position of the g^/ 2 orbital in the HF spectrum, which lies too 



close to the Fermi energy, cf. Fig. 19). The experimental energies of the lowest three states are very close to each 
other, being separated by less than 100 keV, and are well separated by the other two levels. This gross structure is 
already present in the HF result, which, however, greatly underestimates the density of levels. This remains the case 
including pairing correlations at the BCS level (see second column), which lead to a limited improvement. The PVC 
leads to a denser spectrum, considerably improving the agreement with experiment, although one should not expect a 
detailed agreement in the order of three lowest levels. The main effect is already obtained with a single diagonalization 
(NFTO), but the self-consistent treatment leads to an appreciable rearrangement of the spectrum, somewhat reducing 
the initial compression of the levels, slightly improving the agreement with experiment. The increase of level density 
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FIG. 9. The theoretical strength function calculated for the /i u / 2 orbital is compared to the spectroscopic factors associated 
with experimental levels detected in one-neutron transfer reactions. 



can be expressed as an increase of the neutron effective mass from m* rj 0.7m to m* f=a m. One could argue that 
a mean field BCS calculation in a potential with m* — 1 could lead to agreement with experiment in a simpler and 
more direct way. However, such a calculation (still performed with the bare interaction) would greatly overestimate 
the gap [64]. Thus, the simultaneous consideration of gap and low-energy spectra clearly favours a description that 
includes renormalization effects on both quantities. 

The most specific fingerprint of the renormalization processes is represented by the fragmentation of the quasiparticle 
strength, that in our approach is due to the coupling with lqp ® lph states. Experimentally, the effects of such a 
fragmentation leads to a breaking of the observed quasiparticle strength into several peaks, and consequently to 
a reduction of the strength observed in the main peaks, as compared to the full strength expected in mean field 
calculations. In the present case, strength functions measured in one-neutron transfer experiments involving 119 Sn 
or 121 Sn [65-67; show a strong, isolated peak for the quantum numbers J 71 ' = l/2 + , 3/2 + , 7/2 + , 11/2 - . The case 
of 11/2 - is shown in Fig. M where we compare the theoretical values of v 2 (u 2 ) with the observed peak in 119 Sn 
( 121 Sn). Here, and in similar figures, in order to compare the experimental to the theoretical strength functions, we 
have added the energy of the lowest calculated quasiparticle to the experimental excitation energy in the odd nuclei 
119 Sn and 121 Sn. The summed strength of the two main peaks is equal to 0.79 in experiment and to 0.8 in theory. 
The remaining theoretical strength is calculated to lie at higher energy in 121 Sn. A good agreement between theory 
and experiment is found also in the case of the other three quantum numbers, as was shown in Fig. [7] 

It is interesting to notice that the main contribution to the renormalization effects originates from the coupling 
to the lowest vibrational state of each multipolarity : this demonstrates the key role played by the interweaving of 
collective and single-particle modes which is at the basis of the NFT treatment of the elementary modes of nuclear 
excitation. 



This is clearly seen in Fig. 11 



where we compare the full calculations with results obtained including 
only couplings to the lowest vibrational states. It is seen that the coupling to the 2 + low-lying mode is the dominant 
one, providing half of the increase of the gap due to the renormalization effects. The other low-lying modes provide 
about one-third of the remaining increase of the gap, while higher energy modes (mainly giant resonances) complete 
the picture. The strength of the various PVC couplings is reported in Table IhJ where for each pair of orbitals (a, b) 
we list the value of the squares of the PVC matrix elements summed over the phonons of a given multipolarity, 
*^2i V h 2 {ab\v)/(2ji) + 1). We observe that although the deformation parameter associated with the lowest 3~ state is 
the largest one (cf. Tabic m) its influence is hindered with respect to the 2 + by its higher energy and by the fact that 
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FIG. 10. The theoretical strength function calculated for the d 5 / 2 orbital is compared to the spectroscopic factors associated 
with experimental levels detected in one-neutron transfer reactions. 



it acts efficiently only between the d 5 / 2 and the hn/ 2 orbitals, which are the only ones having opposite parity and the 
same spin alignment. We also note that the dominance of the coupling to low- lying vibrational states found in the 
present calculation makes our approach distinctly different from other approaches based on microscopic calculations 
of the PVC with zero-range forces, which need an energy cutoff to avoid an ultraviolet divergence in the self-energy 
[68]. 

The strength of the remaining orbital, 5/2 + , is much more fragmented both in theory and experiment, as is shown 
in Fig. 10 The experiment shows a number of low-energy peaks, which exhaust about 63% of the strength. The 
theoretical distribution is more fragmented than observed in the data; moreover most of the strength lies about 1 
MeV above the experimental findings. In the calculation, the energy of the d 5 / 2 quasiparticle is lowered by its strong 
coupling with si/ 2 state (through the 2 + vibrations) and with the hxxii state (through the 3~ vibrations), as can be 



) 27" and d 



seen in Table fill These couplings bring the energy of the d 5 / 2 state close to the energies of the Si/ 2 

configurations. The matrix elements of the induced interaction Vi n d (cf. Eq. (301) are listed in Table III They are 



3/2 ( 



2+ 



calculated using the values of the renormalized quasiparticle energies E a (l) of the lowest fragment for each orbital. 
Their values can be compared to the typical value of the matrix elements of the bare interaction, which is about 
Go = —0.22 MeV (cf. below Section III B I. They are rather well correlated with the PVC matrix elements reported 
in Table ITTl with the remarkable exception of d§/ 2 : in the latter case, the induced interaction with the orbitals s\/ 2 
and c?3/2 takes large values, associated with the almost degeneracy of the energy of the d§/ 2 quasiparticle with the 
configurations Si/ 2 <g> 2f and d 3 / 2 ® 2+, which leads to the fragmentation of the d 5 / 2 strength. 

These results may indicate that the unperturbed quasiparticle energy of the d 5 / 2 orbital lies too high. In fact, 
moving the energy of the single-particle energy by 600 keV towards the Fermi energy, and leaving everything else 
unchanged, and then solving the BCS and Nambu-Gor'kov equations, one obtains a better agreement with experiment, 
as will be shown below in Section IV A (cf. Fig. [L8J^d)). This shows that a detailed study of the fragmentation process 
can give important indications about how to improve the mean field (cf. also Section VIA). 
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rfg/2 


57/2 


Sl/2 


d-3/2 


fr-ll/2 


^5/2 


0.190 


0.016 


0.467 


0.071 


0.340 


0.088 


0.041 


0. 


0.266 


0.138 


57/2 


0.016 


0.130 


0. 


0.236 


0.020 


0.041 


0.071 


0.168 


0.081 


0.043 


Sl/2 


0.467 


0. 


0. 


0.461 


0. 


0. 


0.168 


0. 


0. 


0.389 


^3/2 


0.071 


0.236 


0.461 


0.253 


0. 


0.266 


0.081 


0. 


0. 


0.092 


^11/2 


0.340 


0.020 





0. 


0.169 


0.138 


0.043 


0.389 


0.092 


0.100 



TABLE II. For each pair (a, b) of the five valence orbitals, we list the value of the sum ^ h 2 (abXv)/(2jt + 1), in MeV 2 . For 
pairs of orbitals of the same parity, the numbers in the upper (lower) row give the contribution associated with A = 2 + (A = 4 + ) 
phonons; for pairs of orbital of different parity, the numbers in the upper (lower) row give the contribution associated with 
A = 3~ (A = 5~) phonons. 





d 5 /2 


57/2 


Sl/2 


rf 3/2 


hu/2 


4/2 


-0.223 


-0.031 


-1.701 


-1.230 


-0.309 


-0.054 


-0.030 


0. 


-0.245 


-0.075 


57/2 


-0.015 


-0.157 


0. 


-0.513 


-0.016 


-0.023 


-0.045 


-0.117 


-0.062 


-0.021 


Sl/2 


-0.383 


0. 


0. 


-0.686 


0. 


0. 


-0.100 


0. 


0. 


-0.177 


d-i/2 


-0.055 


-0.219 


-0.503 


-0.319 


0. 


-0.139 


-0.047 


0. 


0. 


-0.041 


hu/2 


-0.183 


-0.012 





0. 


-0.212 


-0.051 


-0.017 


-0.167 


-0.042 


-0.067 



TABLE III. Matrix elements of the induced interaction Vi n d( a (i)b(i)) (cf. Eq. |30l) between the lowest fragments of the five 
valence orbitals, in MeV. For pairs of orbitals for the same parity, the numbers m the upper (lower) row give the contribution 
of the matrix elements involving 2 + (4 + ) phonons; for pairs of orbital of different parity, the numbers in the upper (lower) row 
give the contribution of the matrix elements involving 3~ (5~) phonons. 
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FIG. 11. The pairing gaps A BCS obtained in the BCS calculation (dots) and A including the renormalization effects (triangles) 
are compared to those obtained including only the renormalization due to the coupling to the lowest 2 + , 3~, 4 + and 5~ QRPA 
modes (filled squares), or only the lowest 2 + mode (open squares). 
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FIG. 12. The BCS (dots) and renormalized pairing gaps associated with the lowest quasiparticle peaks obtained solving the 
Nambu-Gor'kov equations (triangles, cf. Fig. p| are compared to the gaps obtained including the schematic spin-induced 
interaction discussed in the text (diamonds), and to their decomposition into A bare and A pho . 



3. Effects of spin modes 



It is interesting to speculate about the effect of S — 1 modes, that should be included in a more complete calculation. 
Their contribution represents the leading renormalization effect in calculations in uniform neutron matter, where it 
induces a repulsive contribution on the pairing interaction that wins over the attractive contribution associated with 
density modes - thanks to the larger multiplicity of spin modes - leading to a quenching of the pairing gap [6T)]-|72|. 
On the other hand, one expects that the relative weight of spin modes should be considerably less important in finite 
nuclei, where there is no spin equivalent to low-lying surface collective vibrations. 

The repulsive character of spin modes, as opposed to density modes, arises from their different transformation 
properties under time reversal. This is reflected in a sign change in the basic V, W vertices (cf. Eq. (6-207) in [T5]). 
which become (cf. Eq. (18)) 



V(ab(m)\v) 
W{ab{m)\v) 



hs= 
hs- 



'abXv)(u a Ub,m + 
'ab\v)(-u a v btm 



V a Vb, m ) 
f V a Ub,m)- 



(55) 



As a consequence, the product VW which determines the abnormal self-energy Y} 2pho (cf. Eq. (21 )) also changes its 
sign, leading to a repulsive contribution. On the other hand the normal self-energies depend on the square of V and 
W so that their value is increased by the action of spin modes. There is furthermore a change in the expression of 
the basic vertex h, which is given by 



h s =i(abJv) = -(-l) j «- jb l3j u < a\fj(r)\b >< Jb \\(Y x a)j\\ Ja > 



1 



(2j + l)(2J+l) 



1/2 



(56) 



where fj(r) is the form factor associated to the S=l modes transition density. 

The contribution of spin modes to pairing correlations in 120 Sn was estimated in ref. [73], based on a QRPA 
calculation performed with the SkM* interaction. There it was found that, performing a calculation of the pairing 
gap including only the induced interaction, using the simplified expression discussed in Section III A 4 (cf. Eq. (60)), 
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the spin modes decreased the gap arising from the coupling with S = modes by about 25%. We have incorporated 
this effect in our present calculation by introducing a schematic S = 1 response function, consisting of a single peak 
at an energy of about 1Tiuj osc , in keeping with the fact that there are no low-lying collective S = 1 modes in finite 
nuclei. We have used an average value for the product j3j u < a\fj(r)\b >, treating it as an adjustable parameter in 
order to reproduce the 25% reduction in ^, 12 P ho mentioned above, when no bare interaction is included. In this way 
we can estimate the effects of S = 1 modes on the total gap and on the quasiparticle energies and strengths calculated 
in the present paper. 



12 



and 



13 



In 



The values obtained for A and for the quasiparticle strength N are shown respectively in Figs 
keeping with their contribution to the normal self-energy, the action of spin modes tends to reduce, although slightly, 
the value of Z and N, not modifying significantly the agreement with experiment. Spin modes also lead to a reduction 
of A pho by about 25% (150 keV), which reflects directly the reduction introduced in Y} 2pho , while the value of A bare 
is practically unaffected. Altogether, this brings the value of A closer to the experimental value of A exp — 1.4 MeV, 
without deteriorating the agreement with the experimental quasiparticle strength. 
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FIG. 13. The strength iV associated with the main quasi particle peaks in the bare HF+BCS calculation (dots) and solving 
the Nambu-Gor'kov equations (triangles) is compared to the result obtained including the schematic spin-induced interaction 
discussed in the text (diamonds). 
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4- One- quasiparticle approximation and the gap equation 

We have seen that the quasiparticle strength function for the states close to the Fermi energy displays a limited 
amount of fragmentation. If one is only interested in the properties of the main quasiparticle peaks, it may be 
interesting to solve the Nambu-Gor'kov problem by keeping only the main quasiparticle for each orbital in the 
iteration process. One then gets an excellent agreement with the complete calculation for the quasiparticle gap 



A^ p - = Z a 'S a 2,q - p -, as shown in Fig. 14 the result cannot be distinguished from the calculation for A shown in Fig. 

E 



In order to have more insight on these results, we recall the gap equation Eq.(43): 



fi&j = Af cs - £ V ind (anbm)N b{m) S _""> (57) 



Km) y(e b - eF + E— )2 + (s i2 m)) 

While the implementation of the quasiparticle approximation is straightforward, it is important to notice that in 
order to get good agreement one must include the factors N and Z associated with the quasiparticle peaks. Neglecting 
them, that is, attributing the full quasiparticle strength to the single peak retained in the calculation, leads to an 
increase of the gap by about 50%, in keeping with the typical values Z « 0.7, N ~ 0.8. Neglecting N produces 



a smaller effect than neglecting Z, because N affects only the induced part of the gap in Eq. (57). The error is 
especially pronounced for the case of dg/2, where the one-quasiparticle approximation selects the lowest peak, which 
is associated with a small value of N and Z (cf. Fig. Wy. From this result, one concludes that the gap receives only 
very small contributions from the fragments with m > 1, which, on the other hand, account for more than 20% of the 
quasiparticle strength. This is mostly due to the energy dependence of Vi nd , which is strongly peaked at the Fermi 
energy [IT] . 



A gap equation equivalent to Eq. (57) has been obtained in the post scheme (cf. Eq. (54)): 

Aft(m) 



A a(n) = -Z a{n) ]T ^//(o(n)6(m))JV 6(m) ^^-. (58) 

b{m) 



2-E&(rn) 



This expression lends itself to approximations, which are useful in order to make a connection with previous works 
about renormalization effects on pairing. The most important feature of these approximations is the use of simplified 



expressions for the induced interaction Vi nd (cf. Eq. (30)), 

V ind (a(n)b(m)) = J2\, v 2h (2^+i) ) x 

(59) 



E a (n) — Ei 3 (m) — hu)\„ E a (n)+E b (m)+huj\ 

based on the Feynman diagram (lc) in Fig. [I] and neglecting fragmentation: 

j^ (2j b + l) E Q -E a -E b -Tux> Xv 

in which E a = |e — £f\ and Eq is the pairing energy per Cooper pair (equal to about - 2A) [7]. This expression was 
inserted in the BCS gap equation, without taking explicitly into account the renormalization effects on the single- 
particle density, namely, setting the N- and Z— factors equal to one, and using a mean field potential characterized 
by an effective mass m* = m. One then found gaps about 20% larger than those obtained solving the Nambu-Gor'kov 
equations [E]j[H]- The simplified expression for Vi n d was also employed in a more elaborate calculation |12j which used 
HF level from the SLy4 interaction and a constant value N — 0.7. However, £12 was identified with the pairing gap, 
which then tended to be overestimated by a factor l/Z. 

B. Comparison of prior and post scheme with the monopole pairing force 

We have already noticed that solving the Nambu-Gor'kov equations in the post scheme is numerically very demand- 
ing in the case of realistic nucleon-nucleon interactions. In this Section we want to compare the prior and the post 
schemes using a simple monopole pairing force - GqP^ P acting only in our valence space. The present calculation 
is similar to that presented in rcf. [8 , except for details in the mean field (here obtained from a SLy4 interaction 
instead of a Woods Saxon potential with effective mass m*— 0.7), in the QRPA spectrum and in the determination 
of the PVC coupling. We shall however present new results which will extend our previous investigation. 
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FIG. 14. The pairing gaps obtained in the quasiparticle calculation (triangles) are compared to the results obtained neglecting 
the normalization factors N, and neglecting both N and Z. 
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FIG. 15. Comparison of the renormalized pairing gaps A associated with the lowest energy fragments, obtained solving the 
Nambu-Gor'kov equations in the prior (solid line with triangles) and in the post (dashed line) scheme. We also show the bare 
(A 6are ) and the phonon (A pho ) components of the renormalized gaps in the prior (solid lines with diamonds and squares) and 
in the post (dashed lines) scheme. The pairing interaction is a monopole force with coupling constant Go = 0.22 MeV, acting 
in the valence shell around the Fermi energy in 120 Sn, which produces the gap A (solid lines with dots) obtained in the 

BCS calculation. The symbols refer to the position of the various valence orbitals in the SLy4 HF potential. The value of the 
gap A e:E p obtained from the experimental odd-even mass difference is also indicated. 



In Fig. 15 we compare the results obtained in the prior and in the post schemes, for a value Go = 0.22 MeV, 
corresponding to a bare pairing gap jr; 126ore = 1.08 MeV. We have chosen this particular value of Go in order to 
reproduce the average value of the gap A BCS obtained previously with the t>i4 interaction (cf. Fig. rc3j). We see that 
the renormalized pairing gaps are very similar in the prior and post schemes. Moreover the results are also similar to 



27 



Z 

si 



]' 


T 




1 


1 




1 


A-iN 




A 




08 








_ 


0.6 








— 


0.4 


- 


5/2 + 


/ 7/2 + 


l/2 + 3/2 + 




11/2" 


0.2 


- 










- 









i 








2 




-10 




-8 




-6 



e (MeV) 



FIG. 16. Z— and N— factors associated with the various valence orbitals, for the calculation presented in Fig. 
post values are indistinguishable at the scale of the figure. 



FL5l Prior and 



those obtained with the Argonne interaction, shown in Fig. raj not only for the total gap A but also concerning the 
decomposition into bare and phonon contribution. The similarity extends also to the Z— and N— factors, shown in 
Fig. [l6] (cf. Fig. [7k. The fact that the results turn out to be close to those obtained using the Argonne interaction 
is easy to understand in the prior scheme, which is based on the quasiparticle and occupation amplitudes obtained 
with the bare interaction. In fact, the realistic nucleon-nucleon interaction displays a limited state dependence over 
the valence orbitals, with an average value A BCS = 1.1 MeV, very close to the constant BCS gap obtained with the 
monopole interaction with G = 0.22 MeV. 



■ 

3 




i 


i 
• 

■ 


* BCS 

A 

A, post 








- 


— T- 8 
• 








▲ 


A, prior 




M 


• 


• 


- 


>2 




1 


1 


<1 








■ 


• 








1 

r 




11 


■ 
▲ 


* • 

• 








- 


°! 


► 


• 

i 


f 


i 








— J-I 







0.1 



0.2 

G (MeV) 



0.3 



0.4 



FIG. 17. Renormalized gaps A obtained solving the Nambu Gor'kov equations in the prior and in the post scheme with the 
monopole pairing force as a function of the pairing constant Go, averaged over the five valence orbitals. Also shown is the gap 
A scs obtained solving the BCS equation. 
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In Fig. 17 we show the value of the gap, averaged over the five valence orbitals, as a function of the pairing strength 



Go. In the prior scheme, one starts from the u a , v a amplitudes obtained in a previous BCS calculation with the bare 
force, that for G smaller than 0.1 MeV, produces no superfluid solution. On the other hand, the induced interaction 
Vi n d is able to produce a pairing gap by itself. As a consequence, in the prior scheme only £H on contributes to the 
gap, which is therefore underestimated and independent of Go for Go < 0.1. Instead, using the post scheme, the gap 
grows as a function of Go since the bare interaction can provide a contribution to pairing correlations even for small 
values of Go when it is added to the induced interaction through the effective interaction V e f t = V(, ar e + Vind- For 
values of Go larger than 0.1 MeV, the gaps in the prior and post scheme become closer and closer, until the prior 
scheme becomes satisfactory for values of the order of Go = 0.2 MeV. It is interesting to notice that the calculations 
for the bare and renormalized gaps performed with Go = 0.25 MeV reproduce quite well the results obtained using the 
Viowk interaction as bare pairing force, as discussed in Section [VlB| This indicates that Fig. [17] can be used to assess 
in a simple way the effect of renormalization processes on pairing correlations, for the adopted QRPA spectrum. 

IV. CONCLUSIONS 

We have presented a convenient formalism to deal with the basic renormalization processes induced by the coupling 
between quasiparticle and collecive vibrations in superfluid spherical nuclei. We have solved the Nambu-Gor'kov 
equations determining the normal and abnormal energy-dependent self-energies self-consistently. This allows a detailed 
calculation of the low-energy part of the nuclear spectrum in odd nuclei taking into account the fragmentation of the 
quasiparticle strength, as well as the calculation of the pairing gap of the system including the pairing interaction 
induced by the exchange of collective modes. The mean field is based on a Hartree-Fock calculation with the effective 
SLy4 interaction, while the coupling between quasiparticles and vibrations is determined by a QRPA calculation 
that reproduces the empirical polarizability of the low-lying vibrational modes. In the solution of Nambu-Gor'kov 
equations there are then no free parameters, when we use the bare vu Argonne nucleon-nucleon potential as bare 



pairing force. We have obtained a gap equation (cf. Eq. (43)) which takes into account renormalization effects in a 



compact way, and allows one to make contact with previous studies. 

We have discussed in detail the practical implementation of the formalism, introducing two calculational schemes, 
depending on whether the bare pairing interaction and the renormalization effects are taken into account simultane- 
ously (post scheme) or in two steps (prior scheme), performing first a BCS calculation with the bare force. We have 
shown that the two schemes lead to similar results for states close to Fermi energy in the case of a simple monopole 
pairing force with a realistic coupling strength. The prior scheme is computationally much simpler in the case of 
realistic nucleon-nucleon forces which can scatter particles to very high energies. 

We have presented results for the semimagic nucleus 120 Sn, which are in rather good agreement with those extracted 
from one-neutron transfer reactions. The formalism allows also for a clear separation between the contribution of the 
bare force and of renormalization effects: their contributions to the pairing gap turn out to be of similar magnitude, 
confirming previous studies. 

As we have already summarized in Section |IV A| we get an overall agreement with the experimental spectra. This 
agreement appears to be stable, respect to reasonable changes in the ingredients of our calculations (cf . the Appendix) , 
and the contribution of the pairing interaction induced by surface vibration appears to be well established within our 
scheme. 

Several elements remain to be further investigated, before one can reach firm quantitative theoretical results, in 
particular concerning the absolute value of the total pairing gap. On the one hand, the mean field should be either 
derived microscopically, or at least refitted comparing theory and experiment taking into account renormalization 
effects; on the other hand, contributions which are expected to provide a repulsive contribution to the pairing in- 
teraction, associated with three-body effects and with the influence of spin modes, have started only recently to be 
examined. 

A. Hindsight 

One can conclude that it is possible, at least in the case of 120 Sn, to draw a picture that is consistent with the 
available experimental data concerning the quasiparticle properties close to Fermi energy. This is summarized in Fig. 
|18| which shows the results of a calculation with the bare pairing interaction, including all the elements discussed in 
the previous sections. 

The spectroscopic factors and the induced pairing gap are not much affected by the details of the mean field or of 
the effective mass associated with the effective forces we have examined. These quantities are essentially determined 
by the interweaving with low-lying collective surface vibrations, which shifts the quasiparticle energies, increasing the 
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level density close to the Fermi energy, and leads to spectroscopic factors (Z— and N— factors) in the range 0.6-0.8 
(cf. Fig 18 (b)) and to an induced pairing gap A pho of about 0.8 MeV. The coupling of spin modes, due to their 

It reduces the induced pairing 

(a)). 
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higher energy and smaller collectivity, does not affect much the spectroscopic results 

gap by about 0.2 MeV, leading to a final induced pairing gap A pho close to 0.6 MeV (cf. Fig 

The calculated strength function of fragmented levels (the d 5 / 2 in the present case) depends on the detailed position 
of the levels in the mean field potential, and the comparison with experiment provides important information on it. 
In the present case, we have explicitly verified that it is possible to reach a good overall agreement between theoretical 
and experimental spectra, both for the fragmented 5/2 + strength and for the other valence orbitals for which the 
quasiparticle approximation is valid, simply by adjusting the position of the d§/2 orbital in the SLy4 HF potential 
(cf. Figs. 



18 (c),(d)) . 



We also find good agreement between the total pairing gap A and the value derived from the experimental odd-even 
mass difference, starting from the value A BCS = 1.1 MeV obtained with the value of mk associated with the SLy4 



interaction (cf. Fig. 18 (a)). In fact A = A pho + A BCS x Z = 0.6 MeV + 0.8 MeV = 1.4 MeV. 



We remark that microscopic pairing forces lead to a weak state dependence for the orbitals close to the Fermi energy, 
so that in this region we obtain essentially the same results - both at the BCS level and after the renormalization 
process - with a simple monopole interaction, using a suitable value of the pairing constant. Furthermore, different 
values of A BCS , associated with different momentum dependences of mt far from the Fermi energy, essentially lead to 
a shif t of the final gap A, but do not alter significantly the spectroscopic results shown in Fig. 

vTaI 



18 >-d) (cf. Appendix 



Unfortunately, a conclusive theoretical calculation of m^ is not yet available. Furthermore, one should add the 
effects of the three-body force, which is expected to lead to a repulsive contribution to the pairing interaction. 
Interestingly, a recent calculation of A using the V\ ow k pairing interaction and including three-body effects within 
the renormalization group approach, has produced a value for A BCS close to 1.1 MeV in 120 Sn [5DII61J . 
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FIG. 18. Results obtained solving the Nambu-Gor'kov equations including PVC to spin modes with the SLy4 mean field, 
shifting the energy of the ds/2 orbital by 600 keV towards the Fermi energy: (a) renormalized gaps, (b) Z- and iV— factors, 
(c) quasiparticle spectrum, (d) 5/2 + strength function. 
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FIG. 19. Energies of the five HF single-particle levels lying close to the Fermi energy in Sn, calculated with the effective 
interactions SKa, SLy4 and SGII. The effective mass associated with these forces at saturation density is also indicated. 

VI. APPENDIX 

In this Appendix we consider the sensitivity of our results to some of the prescriptions and ingredients of the 
calculations. 



A. Mean field 

The results obtained in the main text have been obtained with a HF mean field produced by the SLy4 interaction, 
whose effective mass at saturation is rrik = 0.7m. In this Section we show results obtained with two different 
interactions of the Skyrme type: the SKa [75] and the SGII [77], having respectively a lower (mu — 0.61m) and 
a larger (m^ = 0.78m) effective mass. All the other features of the calculations are the same as for SLy4. In 
particular, the PVC matrix elements are the same, because in our approach they only depend on the properties of the 
phcnomenological phonons and are calculated with the wavefunctions of a fixed potential with effective mass equal 
to one. As expected, the SGII interaction produces a mean field associated with a higher level density, leading to a 

the value of A BC is equal 
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larger pairing gap in the BCS calculation with the V14, potential, as shown in Fig. 

to about 1.6 MeV, to be compared with the value 1.1 MeV previously obtained with the SLy4 interaction (cf. Fig. 

|6J). The value obtained with the SKa interaction is instead equal to about 1 MeV (cf. Fig. |20|fo)). 

The renormalization processes act in a very similar way as was previously found in the case of SLy4. The average 
value of the phonon induced gap A pho is in all cases equal to about 0.8 MeV: in the case of SGII A pho accounts for 
about 30% of the total gap, while in the case of SKa A pho and A 6are are of similar magnitude, as in the case of SLy4. 
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FIG. 20. The st ate -dependent paring gap A calculated in BCS with the bare V14 interaction is compared to the renormalized 
gap A (cf. Eq. (381 ) obtained solving the Nambu-Gor'kov equations. We compare results obtained with a mean field produced 
with the SGII interaction (a) and with the SKa interaction (b) (cf. Fig. [6] for the corresponding calculation with the SLy4 
mean field) [75]. The symbols refer to the position of the various valence orbitals in the SGII HF potential. We also show the 
decomposition of A into the bare and phonon contributions A 6are and A pho . The value of the Fermi energy tF and of the gap 
obtained from the experimental odd-even mass difference A exp are also indicated. 



The low-lying spectra are shown in Fig. 21 while the Z— and N— factors of the lowest fragments are reported 
in Fig. [23] Also in this case, the action of the PVC is similar to that already discussed in the case of the SLy4 
interaction (cf. Fig. pi. However, the agreement with experiment is better than obtained previously. This is due to 
the initial position of the d 5 / 2 and g^/2 orbitals, which are closer to each other; furthermore, the d 5 / 2 single-particle 
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FIG. 21. The theoretical quasiparticle spectra obtained at the various steps of the calculation are compared to the experimental 
data. We compare results obtained with a mean field produced with the SGII interaction (a) and with the SKa interaction (b) 
(cf. Fig. [8] for the corresponding calculation with the SLy4 mean field). 



calculated with SGII lies at a distance of about 2.2 MeV from the Fermi energy, instead of 2.8 MeV as in the case of 
SLy4: this leads to a good description of the fragmentation of this orbital, as shown in Fig. [22] that in the case of 
SLy4 was obtained only shifting the energy of this level (cf. Figs. 10 and 18 1. The quality of the spectrum obtained 
with SKa is almost as good as with SGII: however, the orbital <?7/ 2 lies more distant from the Fermi energy and 
becomes fragmented, contrary to experiment (cf. Fig. [23]). It is interesting to notice that taking the SLy4 mean field, 



and changing the energies of the five valence orbitals with those associated with the SGII interaction, one obtains 
practically the same spectrum and the same Z— and N— factors as with SGII. This in spite of the fact that A BCS 
becomes equal to about 1.25 MeV, to be compared with 1.1 MeV (SLy4) and 1.6 MeV (SGII). At the same time, the 
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average value of A is equal to about 1.7 MeV, to be compared with 1.6 MeV (SLy4) and 2.1 MeV (SGII) (cf. Fig. 



24 1. These results show that the low-energy spectrum is determined by the position of the valence orbitals, while the 
absolute value of the BCS gap also depends on the effective mass associated with distant levels. 

We can conclude that rcnormalization effects arc similar for the three mean fields we have considered. The compar- 
ison with the odd-even mass difference favours forces having low effective mass (SLy4 or SKa). The PVC improves 
the agreement of the spectral properties with experiment. However, the quality of this agreement depends on the 
specific position of the mean field single-particle levels close to the Fermi energy. The value of the effective mass far 
from the Fermi energy determines the magnitude of the final gap A by shifting the value of A BCS , while the value of 
A pho and the properties of the low-lying spectrum are not very sensitive to it. 
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FIG. 22. The theoretical strength function calculated for the d 5 / 2 orbital is compared to the spectroscopic factors associated 
with experimental levels detected in one-neutron transfer reactions. We compare results obtained with a mean field produced 
with the SGII interaction (a) and with the SKa interaction (b) (cf. Fig. 10 for the corresponding calculation with the SLy4 
mean field). 
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FIG. 23. Comparison of the N— and Z— factors associated with the lowest quasiparticle peaks in the Nambu-Gor'kov 
calculation shown in Fig. 20 We compare results obtained with a mean field produced with the SGII interaction (a) and with 
the SKa interaction (b) (cf. Fig. FT] for the corresponden calculation with the SLy4 mean field). Also shown by stars are the 
values of the experimental quasiparticle strength 66], except for the d$/2 orbital which shows a pronounced fragmentation. 
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FIG. 24. The same as in Fig. |20| for the pairing gaps calculated with the SGII mean field, but substituting the single-particle 
energies of the five valence orbitals with the energies calculated in the SLy4 mean field. 
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B. Bare pairing interaction 



The calculations reported in the main text have been carried out adopting vn as the bare pairing force. However, 
the Viowfc version of the Argonne potential has been used by several groups and here we show the results using this 
bare pairing force in our Nambu-Gor'kov formalism. The corresponding bare and renormalized gap are shown in Fig. 



25] The average value of the bare gap is equal to about 1.4 MeV, in agreement with ref. [57], to be compared with 
the value 1.1 MeV obtained with the Argonne interaction. The effect of the renormalization processes increases the 
gap on average by about 500 keV, similar to the case of v±4, (600 keV, cf. Fig. [6]). Also the values of the phonon 
induced component A pho are quite similar. The results of the calculations with Vn and Vi ow k turn out to be quite 



similar to those obtained in Section III using the monopole force respectively with coupling constants Go = 0.22 and 
G = 0.25 MeV (cf. Fig. fTf}. 




£ (MeV) 



FIG. 25. The state-dependent paring gap A ca lcul ated using the Viowfc potential with a cutoff A = 4 f m x as pairing force 

is compared to the renormalized gap A (cf. Eq. (38 1 ) obtained solving the Nambu-Gor'kov equations. The symbols refer 
to the position of the various valence orbitals in the SLy4 HF potential. We also show the decomposition of A into the bare 
and phonon contributions A 6are and A ph °. We thank A. Pastore for providing us with the HFB calculation with the V\cmk 
potential. 
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C. QRPA 



We recall (cf. Section [n]) that the results shown in the main text are based on effective PVC vertices calculated with 
single-particle levels and wavefunction associated with a Woods-Saxon potential, with an associated effective mass 
m* — m. Another reasonable choice, often adopted in the literature, would be to use instead the same single-particle 
levels used in the HF+BCS calculation. Using this prescription and readjusting the coupling constants \\ °f the 
QRPA calculation so as to reproduce the experimental properties of the low-lying vibrational states, we obtain the 
renormalized pairing gaps shown in Fig. |26| It is seen that the difference is not critical, our preferred choice leading 
to gaps which are smaller by about 10%. 



2.5 



> 

1 


1.5 


< 


1 




0.5 









12 



5/2 



7/2 + l/2 + 3/2 + 



11/2 



BCS 



A (SLy4QRPA) 
A (WS QRPA) 



10 -8 

8 (MeV) 



-6 



FIG. 26. The renormalized pairing gaps A, calculated using the levels of a Woods-Saxon potential to compute the QRPA 
spectrum and the PVC vertices, are shown by triangles (cf. Fig. pj and are compared to the results obtained using SLy4 
single-particle levels (squares). We also show by dots the pairing gap A scs calculated in BCS. 
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FIG. 27. While in the rest of the work we have considered the effect of the coupling betwene quaisparticle and collective QRPA 
phonons (left), in this section we shall consider the coupling to uncorrelated excitations (right). 
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FIG. 28. Pairing gaps associated with the main quasiparticle peaks obtained in the one-bubble approximation. We show the 
total renormalized gap A and its decomposition into the bare and the phonon contributions. In the inset, the Z— values are 
also shown. 
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D. Bubble approximation 

It is interesting to assess the importance of the collectivity of QRPA modes in the pairing calculation. This can 

be. doiiebYXompajdng.ouiTffl-evious^rcswife. with calculations .performed coupling the quasiparticles .with unperturbed 
p-h excitations \ci. rig. 27). In Pig. 28 we show the pairing gap and tire Z— factors calculated at the p-h level: 



they are very close to 1, indicating that the renormalization processes are much less important as compared to the 
calculations with collective modes. Correspondingly, the phonon contribution to the gaps is drastically diminished 
(compare with Figs. [6] and [7]). Also the quasiparticle spectrum becomes very close to that obtained at the HF+BCS 
level. On the other hand, the total gap is reduced by only 10% as compared to the QRPA result: this happens 
because, even if the phonon-induced interaction is considerably weaker, the product A = ZA BCS + A pho is much 
less affected, because the Z— factor is correspondingly increased. Schematically, while in the QRPA case one finds 
Z w 0.75 and A = 0.75 x 1.1 + 0.8 = 1.6 MeV, in the present case Z « 0.95 and A = 0.95 x 1.1 + 0.3 = 1.4 MeV. 
One then concludes that, although the value of the gap is similar using the QRPA or the unperturbed response, the 
physical picture is quite different. This is related to the fact that the contribution of the phonon interaction to the 
total gap is equal to about 50% considering the coupling to QRPA modes, and to about 20% coupling to unperturbed 
p-h excitations. Furthermore, we have performed a calculation correcting for the Pauli principle, namely removing 
the contribution from a bubble when the particle or the hole coincides with the intermediate state {nlj}, with the 
proper angular momentum factor l/(2j + 1) (cf. [TJj, Appendix F). The pairing gap is reduced by less than 5%. This 
represents an upper limit for role of the Pauli correction in the calculations based on collective phonons reported in 
the main text. 
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